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^ r*| Abstract 

Using our recently developed one parameter quasiparticle model, we analyze more recent (refined) 

O l' results of (2+1) flavor QGP in lattice simulation of QCD by various groups [T] [2] [3] [4]. We got a 

remarkable good fit to lattice thermodynamics of Ref. [I] and reasonable good fit to Ref. [2] by adjusting 

i-Q ■ single parameter of the model which may be related QCD scale parameter. Further we extend our model 

for above system with zero chemical potential to non-zero chemical potential and predict quark density 

without any new parameters which may be compared with future lattice data. 
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1 Introduction : 

Recently, there are various attempts to study lattice quantum chromodynamics (QCD) with quarks by 
simulations at finite temperature and chemical potential which describes a system with quarks and gluons in 
plasma state, called quark gluon plasma (QGP), at thermodynamic equilibrium and various thermodynamic 
quantities are derived. Thermodynamics of QGP is important to describe the evolution of fire-ball formed 
in relativistic heavy ion collisions, astrophysical objects and early universe so on. It is a challenging job 
to minimize lattice artefacts and obtain a reliable results. Eventhough results of various groups [TJ ® [3J 0] 
predict nonideal effects qualitatively, there were differences in quantitative results because of lattice artefacts. 
Therefore, it may be useful to analyze lattice gauge theory (LGT) results of various groups with theoretical 
models. 

Quasiparticle model is one of the simplest model, widely studied in literature [6j [7] . There are varieties 
of quasiparticle models, but all of them with few phenomenological parameters. We have also developed 
a quasiparticle model with only one parameter [5]. We use this model to study, in a unified manner, the 
LGT results of various groups and a reasonable good fits were obtained to all results by adjusting the single 
parameter which may be related QCD scale parameter. Difference between various lattice results reflects as 
different values for the parameter of the model. 

Next section we discuss the concept of statistical mechanics of quasiparticle model. In section III quasi- 
particle model for (2+1) flavor QGP is discussed. Results and conclusions follows in section IV. 

2 Quasiparticle Model: 

It is a phenomenological model to describe the statistical mechanics and thermodynamics of a non-ideal 
system. It is assumed that the effects of interaction leads to thermal excitation of quasiparticles with thermal 
mass which is related to interaction strength or coupling constant and temperature. We assume that whole 
thermal energy is used to excite this quasiparticles and it is an ideal system of quasiparticles. It means 
that system with real particles with mutual interaction at thermodynamic equilibrium is approximated by 
ideal system with quasiparticles at thermodynamic equilibrium. However, all results of standard statistical 
mechanics of ideal system may not be applicable here as such, because of the temperature dependent thermal 
masses of quasiparticles. 



Suppose we consider such a system in canonical ensemble, then we know that the system is in thermal 
contact with a reservoir and attains the temperature of reservoir at equilibrium. Partition function is [S] 

r 

and average energy, by definition, 

U=<E r >=-^, (2) 

where (3 = 1/T and E r is the r th energy state of the system. In quasiparticle model E r also depends on 
T through the mass of quasiparticles and Eq. @ is just an identity, namely, explicit derivative of partition 
function with respect to /3. In terms of T, above equation becomes 



dT Qn ^ dT 

which may be written as 



T2 ^X^ Er % ^ 



tj _ j,2 d l n QN ^2 ] V^.-SR^t 



U d ( „ f T , 2 dm( T ) 8E r \ 



Other thermodynamic quantities are obtained from the Helmholtz free energy, 



BA 
A = U-TS = U + T—, (5) 



which may be written as 

U _ d A, 

T 1 ~ dT^T' 

Comparing with Eq.@, we get, 



(6) 



A = -TlnQ N+ [ T drr^< d ^>. (7) 

Jt ot dm 

This is known as the bridging relationship between statistical mechanics (right hand side) and thermody- 
namics (left hand side). Since the Helmholtz free energy is a thermodynamic quantity, it is a function of 
temperature (T), Volume (V) and number of particles (N) and supposed to be unaware of the temperature 
dependence of quasiparticle mass. Quasiparticle mass enters only through statistical mechanics, right hand 
side of Eq.([7]). All other thermodynamic quantities may be obtained from A(T, V, N) using thermodynamic 
relations. For example, pressure is given by 

dA 

P=- W =P l d + f(T), (8) 



where pid is the expression for pressure of ideal system. Thus we get an extra temperature dependent term 
f(T) to pressure compared to ideal gas. Many earlier authors [B] and even today [TU], overlook this term and 
neglect f(T) and lead to thermodynamic inconsistency and so on. Of course, there are attempts by [llj to 
unify all different models in a thermodynamically self-consistent models. We still believe that our approach is 
a natural extention of normal statistical mechanics to temperature dependent quasiparticle system, mainly 
because left hand side of Eq. (|4l [7]), being a thermodynamic function, should not explicitly depend on 
m(T), but over all dependence on temperature. Whereas statistical mechanics, right hand side depends on 
m(T). Therefore, in our analysis we first evaluate U or e = ¥?, energy density and derive pressure and other 
thermodynamic quantities using thermodynamic relations. 

3 (2+1) flavor QGP: 

This system is assumed to be consists of gluons, up and down quarks, and strange quarks where all of them 
treated as non-interacting quasiparticles with respective temperature dependent thermal masses in addition 
to physical masses. Following the standard text book statistical mechanics [9], we calculate e(T) first and 
then p(T) is obtained from thermodynamic relations. Statistical mechanics and hence m(T) enter only in 
e(T). Treating as grand canonical ensemble with zero chemical potential we have, 



T4 2 Wo gV^f-l' [ ' 



for gluons and 



d roc x 2 Jx 2 +m 2 
% = £ dx V \ (10) 

T 2n Jo e v x + o + 1 



for quarks. d g = 16 and d q = 12 are the internal degrees of freedom of gluons and quarks. rh g = m g /T and 
fh q = rriq/T, where m g and m q are the total masses of gluons and quarks, respectively, including thermal 
masses. For (2+1) flavor QGP, total energy density is 

e = e g + 2e u + e s , (11) 

where u and s refers light (up and down) quarks and strange quark, respectively. Gluon mass may be taken 
as 01], 

m 2 g = \<4, (12) 



and quark mass 

m 2 q = (m q0 + mf) 2 + mj (13) 

Plasma frequencies u p and to./ arise because of the collective behavior of plasma due to mutual interactions 

among particles. In the absence of interaction, whole thermal energy may be used to excite real particles 

from the vacuum. From the perturbative QCD calculations we have 

2 T 2 2 T 2 

w 2 = £_(2JV c ) + iL_n / . (14) 

where N c = 3 for QCD and n/ is the number of flavors, g is the coupling constant. Similarly, 

_ {Nl-l) g*T> 

m f~ 2N C 8 ' (15) 

Above expressions for plasma frequencies are valid at high temperature. To include low temperature as well, 
we model them as, 



2 2 2 "5 V~^ 2 2 ">q na\ 

u p = a g g j7 + }^ a i9 V' ( 16 ) 



and 



2 J.2 2 "? /-, rr\ 

such that at high temperature they goes to QCD perturbative results Eqs. flTll [T5)) and the coefficients a g , 
a q and & g may be determined. n g and n q are densities of gluons and quarks respectively and defined as, 

r h_ = ^9_ f°° dx x ( 18 \ 



■' .---??- / dx - _ , (19) 



for gluons and 

r 3 "27r 2 7o """eV^+^ + l 
for quarks. 

Above form for frequencies is motivated by similar work in ultra- relativistic (e, e + ,7) plasma [121 113 ) . 

2 

where plasma frequencies are proportional to 2 ^- . This method works well to explain LGT results of pure 
gluon, 2 flavor, 3 flavor and earlier (2+1) flavor QGP LGT results [8]. S§ and S| are obtained from it's 
definition Eqs. (|18l I19[) which depends again on thermal masses and hence becomes integral equation which 
need to be solved self-consistently and which on substitution in Eq. (jlip we get energy density. Then all 
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Figure 1: Plot of e/T 4 and P/T 4 as a function of T for (2+1) flavor QGP from our quasiparticle model and 
compared with LGT results of Ref. jj~] (symbols), Fig. 1(a), and Ref. [2], Fig. 1(b). 

other thermodynamics may be obtained from thermodynamic relations. We use 2-loop running coupling 
constant, 



°.(T) = fc 



67T 



3(153 - 19rt/) ln(21n(T/A T )) \ 
4tt (33-2n/)ln(T/A T ) V (33 - 2n/) 2 ln(T/A T ) J ' 



(20) 



which appear in thermal masses. Only parameter in the model is the At which we adjust to get the best fit 
to LGT results. 

4 Results and Conclusions: 



First, we have to solve equations for densities of up quarks, strange quarks and then gluons, separately, which 
are all integral equations, Eq. (|191 [TS]) . and then these densities are used to get the thermal masses, Eq. 
(fT6l [P7f . Using these thermal masses energy density, Eq. ((4]) , is calculated. Only parameter in the model 
is At which we adjust such that we get good fit to LGT results. However, we found that LGT data of [5] 
lies above our results and that of [I] lies below our results. Hence, we multiply LGT data with appropriate 
factor, as done in earlier quasiparticle models models [7 which may be accounted for uncertainties in LGT 
data. In Fig. 1(a), we compare our results with LGT results of Ref. [T|. Physical masses of strange quark is 
taken as 150 MeV and light quark mass is 28 1 15 of strange quark mass. To get good fit we need At = 135 
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Figure 2: Plot of e/T 4 as a function of T for (2+1) flavor QGP from our quasiparticle model and compared 
with LGT results (symbols) of Ref. [3], Fig. 2(a) and Ref. [4], Fig. 2(b). 



MeV and LGT data need to be multiplied by a factor 1.07. Similar comparison is made for LGT data of Ref. 
[5] in Fig. 1(b) and good fit is obtained for At = 80 MeV and with LGT data multiplied by a factor .99. 
Strange quark mass is taken as 150 MeV and light quark mass is ^ of strange quark mass. We see that two 
different LGT data lie in two opposite sides of quasiparticle model. The discrepency in two different LGT 
data was attributed to uncertainties associated with lattice artefacts. For both the LGT results we took the 
respective parameterized functions presented in references [3 [2] which reproduces LGT data. In Fig. 2, we 
also compare our model with preliminary results of [3] and [4], and fitting parameters are At = 150 MeV 
with multiplication factor 1.2 and Ax = 105 MeV without any multiplication factor, respectively. Pressure 
(P/T 4 ) is obtained from energy density using thermodynamic relations 



dP 

e = T- P 

dT 



(21) 



on integration over temperature as done in Ref. [8]. Since quasiparticle model is valid for plasma state, 
lower limit of integration is taken to be critical temperature. We see from Fig. 1(a) that quasiparticle model 
fits LGT data of Ref. [T] surprisingly well for both the energy density and pressure. In Fig. 1(b) we have 
similar comparison with LGT results of Ref. ,2_ and obtain a reasonably good fit except very near to T c . 
For other LGT data our results are compared in Fig. 2(a) for Ref. [3] and Fig. 2(b) for Ref. 0]. In these 



two cases we don't have any parameterized functions to represent LGT data to fine tune the parameter of 
our model. 
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Figure 3: Plot of e/T 4 and P/T A as a function of T for (2+1+1) flavor QGP from our quasiparticle model 
and compared with LGT results of Ref. 1 (symbols) . 



It is interesting to note that our model, extended to (2+1+1) flavor QGP by including charm quarks 
with physical mass equal to 11.85 times the strange quark mass, also fits the corresponding LGT data of 
Ref. |T] reasonably well as shown in Fig. 3. Here At = 135 MeV and multiplication factor is .93. 

So far we concentrated on a system with zero chemical potential and next we extend the model to system 
with finite chemical potential. Following the literature [3 [8], we just need to replace T/ At in coupling 

constant, Eq. (|2U)) . by 

T I //2~ 

(22) 



T 
A T 



— Wl + (1.91/2.91) 2 ^, 



j^2 



using the results of Schneider |14[ [5] and replacing perturbative plasma frequencies by, 



g 2 T* 



(i 



7T 2 T 2 



) , 



and 



^ = 4S(( 6 + n /) + ^#), 



(23) 



(24) 



p 18 u /; 7T 2 T 2 

where fif is the quark chemical potential of different flavors. Here we consider a system with \x u = fi c i = /j 
and all other chemical potentials are zero. With these modifications, but without any new parameters, light 



quark density, N q = (nj — n{), is plotted in Fig. 3(a) and Fig. 3(b) corresponding the fitting parameters 
of Ref. pQ and Ref. [2] respectively which may be compared with lattice results in future. There is a wide 
difference between two results except at very high temperature, m, fi\ refers to light (up and down) quark 
and antiquark densities. 
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Figure 4: Plot of N q /T 3 as a function of T for (2+1) flavor QGP from our quasiparticle model with the 
fitting parameters of Ref. [I], Fig. 3(a) and [2], Fig. 3(b) for fi/T c equal to 1.0, 0.7, 0.4 and .1 (top to 
bottom curves respectively). 



In conclusion, we applied our one parameter quasiparticle model to (2+1) flavor QGP system and studied 
the statistical mechanics and thermodynamics, and compared with lattice data of various groups. It is re- 
markable that just by adjusting one parameter, QCD scale parameter of 2-loop coupling constant, reproduces 
lattice results of various groups and it's values lie between 80 to 150 MeV. No need of use of full LGT data 
of one of the thermodynamic quantity to obtain model function, as done in some quasiparticle models, and 
then evaluating other thermodynamic quantities and fitting to corresponding LGT data. This corresponds to 
infinite parameter model and definitely will fit LGT data because of thermodynamic relations. Shortcomings 
of our model may be the ignorance of the QCD parameter At and the origin of multiplication factors which 
are different for different lattice groups. Of course, almost all other models used to describe the non-ideal 
behaviour of QGP have more than one parameter. 
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